Spin Transport in a Quantum Wire 
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We study the effect of electron-electron backscattering interactions on spin transport in a quantum 
wire. Even if these interactions have no significant effect on charge transport, they strongly influence 
the transport of spin. We use the quantum Boltzmann equation in the collision approximation to 
derive equations of motion for spin current and magnetization, fn the limit of small perturbations 
from equilibrium, we explain the existence of 'precessional' and 'diffusive' behaviors. We also discuss 
the low-temperature non-linear decay of an uniform spin current outside the hydrodynamic regime. 
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I. INTRODUCTION 

Recent developments in the study of one-dimensional 
(ID) physics have been aimed at observing the somewhat 
intriguing processes of spin-charge separation and charge 
fractionalization. The theoretical study of various inter- 
acting ID systems predicts the decay of the electron into 
separate spin and charge excitations traveling at different 
velocities (see e.g. |lj). Also, the charge of the excitations 
is expected to renormalize in the presence of interactions. 
This phenomenon is called charge fractionalization . 
This strongly non-Fermi liquid behavior is described by 
the Luttinger Liquid theory || . 

With the discovery of quasi-one-dimensional systems 
such as quantum Hall effect edge states, quantum wires, 
and more recently carbon nanotubes, there has been good 
experimental confirmation of this theory How- 
ever, direct observation of charge fractionalization or 
spin-charge separation remains a challenge, and it is be- 
lieved that spin transport may provide useful insight into 
the detection of the theoretically predicted separation 
between charge and spin [p|-^0[. For example, as seen 
here, backscattering interactions terms markedly affect 
spin transport, but do not affect charge transport. 

Besides providing insight into the physics of Luttinger 
Liquids, spin transport has also useful technological ap- 
plications Illj]. Because of the long lifetime of spin ex- 
citations, the spin of the electron may be useful for in- 
formation storage or as a transport element, and could 
have significant applications in fields such as quantum 
computation |]l2| . 

This paper investigates the effect of electron-electron 
interactions on spin transport in ID. We consider the 
model of an ideal, infinitely long quantum wire with two 
transport channels: spin up and spin down. As is well- 
known, the forward scattering interactions are responsi- 
ble for charge fractionalization and spin-charge separa- 
tion, and significantly affect charge transport. By con- 
trast, the backward scattering processes are marginally 
irrelevant under the renormalization group flow |13|, and 
do not contribute to the transport of charge. However, 



the effect of backscattering processes on the spin dynam- 
ics can be quite dramatic, ft yields two different regimes 
S: one characteristic for small temperatures where the 
magnetization vector spatially precesses about the direc- 
tion of the spin current, and another one relevant for 
higher temperatures where spin diffusion dominates. In 
our analysis we will only take into account the effects 
of the backscattering processes; due to spin-charge sep- 
aration, we expect these results to be unchanged by the 
forward scattering. 

Equations of motion for the spin current and magne- 
tization have been derived in Q, mostly phenomenolog- 
ically. Here we provide a more general approach. The 
treatment we use follows closely the Kadanoff-Baym for- 
malism Ji4| , generalized to include the electron spin and 
chirality. We use the quantum Boltzmann equation in 
the collision approximation to write down a general form 
for the spin transport equations. 

In the hydrodynamic limit of small perturbations from 
equilibrium, we retrieve the behavior predicted in . In 
this limit we can identify the terms responsible for the 
precession of the magnetization and for the spin diffusion. 
We find the spin diffusion coefficient 
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in qualitative agreement with the estimate of Ref. S , but 
with a different numerical prefactor. This linearized hy- 
drodynamics is manifestly spin-charge separated, giving 
us some confidence in the validity of the result. 

In the presence of larger perturbations, non-linear 
terms appear as well. An analysis of the full equations 
of motion may be done numerically. For the particular 
case of an injected uniform spin current, we perform a 
numerical analysis of the equations at various tempera- 
tures. At high temperatures, we retrieve purely diffusive 
behavior. In the limit of zero temperature, however, the 
system ceases to show exponentially decaying behavior 
and we can find an analytic solution for the decay of the 
spin current. This solution is in very good agreement 
with the corresponding numerical results. 
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We check that the backscattering interactions cause no 
corrections to the charge and heat transport equations. 

A comment on the validity of these results is in or- 
der. It seems quite likely that the quantum Boltzmann 
approach is asymptotically exact in the hydrodynamic 
regime, due to the marginal irrelevance of backscatter- 
ing in equilibrium (up to possible logarithmic corrections 
requiring a higher-order analysis). The low-temperature 
limit is, however, much less transparent. Indeed, the 
quantum Boltzmann description is almost certainly not 
exact in this regime. In particular, while the original for- 
mulation is manifestly spin-charge separated, the quan- 
tum Boltzmann equations are not, but become so in lin- 
ear response. Thus the decay of the spin current found at 
low temperatures depends upon the form of the charge 
distribution functions defined in Section IV. Neverthe- 
less, this dependence appears fairly innocuous, and leads 
to no obvious unphysical behavior for local quantities 
(e.g. magnetization, charge density, current, energy den- 
sity). For this reason, and since the results are in agree- 
ment with expectations based on scaling, we find the 
quantum Boltzmann approximation reasonable an this 
limit as well. 

In the following section we write down the spin trans- 
port equations. The derivation of these equations is pre- 
sented in detail in the Appendix. In section III, we study 
the perturbative regime and the resulting precessional 
and diffusive behaviors. In section IV, using numeri- 
cal and analytical methods, we study the time evolution 
of a uniform spin current injected into the system. We 
present the conclusions in section V. 

II. SPIN TRANSPORT EQUATIONS 

The non-interacting low energy effective Hamiltonian 
for a spin-i quantum wire is given by: 

H ° = X dX hv F^La id ^L a ~ ^rJ^RvI C 1 ) 
a J 

where a =f / J, is the spin index, R/L denotes the chi- 
rality of the right- and left- moving modes, and vf is 
the Fermi velocity. The quantum wire is taken to be 
infinitely long. We consider adding two types of local 
interactions that are allowed by SU(2) invariance. The 
forward scattering interaction is described by: 

Hi=-u s I dxJ R (x)J L (x) (2) 

while the backward scattering interaction is: 

H 2 = -u J dx3 R {x) ■ J L (x) (3) 

= -| J dx ^R a { x )^ocp^Rp{x) ■ ^[ 7 {x)a lS ipL8(x). 



Of the quartic processes described by Hi and H 2 , only 
the backscattering processes will contribute to spin trans- 
port. The forward scattering terms are included as a 
check. 

Given the above form for the interactions, we are in- 
terested in deriving equations for the generalized Wigner 
distribution functions f^ ia2 (p;x,t). For each chirality, 
these functions form a 2 x 2 matrix in the spin space. 
If we choose the z axis as the reference axis, the diago- 
nal components of this matrix give the distribution func- 
tions for the particles with the spin directed upward and 
downward along this axis. The off-diagonal components 
are related to the distributions of particles that have the 
spin directed along the x and y axes. In equilibrium the 
off diagonal components vanish, while the diagonal ones 
are equal and independent of space and time, and the 
Wigner matrix can be written as: 

ff{p V x,t) = f r (p)5 a p, (4) 

f°R/L(p) = —n^ • (5) 

e k B T + 1 

We can relate the Wigner matrices to the actual spin 
current and magnetization, as well as with the total 
charge density and charge current. The local densities 
of the charge and of the spin components in the right 
and left moving channels are given by: 

KW = \ E / J? ia, (p^C ( 6 ) 

Ql,Q 2 •'P 

n r (X) = J2 I fr a (p;X), (7) 

a J P 

and the total spin current and magnetization are 
J s = v F (M R ~M L ), 

M = M R + M L . (8) 

For simplicity we chose to denote (x, t) by X and J |^ 
by J and we will make use of this notation throughout 
the paper. 

The total charge density and current can be written as 
I = v F (n R - n L ), 

p = n R + n L . (9) 

In equilibrium all the components of the spin current 
and magnetization vanish. If the system is perturbed, the 
effect of the interactions is to move the system towards 
its equilibrium position. The equations that describe this 
dynamics can be derived using a technique similar to [ pT[ . 
The details of the derivation are presented in the Ap- 
pendix. 

The resulting generalized Boltzmann equations for the 
Wigner distribution functions /" lQ2 can be written as: 



(d t + VF d x )f^(p)= ^ur(p) J q f2 a2 (i) j /r i7 (?)/r 2 (p)] 



(10) 
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f2 3a3 (i)-fR a3 (p) 



fR 3a3 (p)-fl 3a3 (q) 



while the corresponding equation for /^ lQ2 can be obtained from Eq.(10) by interchanging L and R everywhere and 



by replacing (d t + Vpd x ) by (d t — Vpd x ). To simplify the notation, as all the distribution functions are evaluated at 
the same position x and at the same time t, we dropped the spatial and temporal coordinates. 



III. PERTURB ATIVE EXPANSION integral form of Eq.([Io|) that is given in the Appendix 

in the equation ( |A2fj| ). We expand this form in F" 13 , 
If the deviations from equilibrium are small, we can so that in the terms proportional to u 2 we keep only 
write the matrices of Wigner distribution functions as a the linear terms in F. We make use of the fact that 
sum of equilibrium diagonal matrices and small pertur- Ir{p) + flip) — 1 an d we obtain 
bations, f?P(p;X) = f?(p)6 a0 + F?P(p;X). We use an] 



(d t ±v F d x ) J /$2 a (P;*)=±|ri£ J F^(q;X) J F^ (p; X) T i± £ / ^ *) J ( n ) 
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s, 

Ah vf 



\F^( q - : X)-F^(q;X)] f f R (p)f L (p) 

Jp 

'ia 2 X! / 
J a 



F^(q;X)-F^(q;X) / f R {p)f L {p). 



R ± L 

a 3 J 1 - " J P 



We use Eq.(g) and 

'/£(p)/£(p)= 



dp 



n Rpvp -hpv F 

e k B T -fie k B T 
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(12) 



to write 



(ft ± v F d x )M R/L (X) =±^E M^(*)M* (*)e>" fc 



3>fc 



u 2 /c B T 
~^h 3 v 2 F 



[W R (X)-Ml(X)}. 



(13) 



The temperature dependence of the second term only 
is a result of an explicit dependence of this term on the 
form of the equilibrium distribution functions, while the 
first term linear in u is independent of this form. 

In terms of the total spin current and magnetization 
Eq.(Oh becomes: 



(14) 



d t 3 s (X) + v 2 F 8 x M(X) = -M(X) x J S (X) 

h 

2 k B T T fV\ 

2-kK v F 



We also note that the spin continuity equation still 
holds: 



d t M(X)+d x J s (X) = 0. 



(15) 



Equation (JIJ) is qualitatively similar to the corre- 
sponding result in S, up to a numerical factor in front 
of the J s term. 

The charge dynamics is unaffected. We obtain from 
Eq S .(M): 



d tP {X) + d x I(X) = 0, 



d t I{X) + vpd x p(X) = 0. 



(16) 



(17) 



As expected, this result holds as well in the non- 
perturbative regime, and it can be exactly derived using 
the more general equation (A20). 

We analyze Eq.(|l4|) to extract the physics of the sys- 
tem. If the first term of the RHS of the equation domi- 
nates, as it is the case for zero or small temperatures, the 
spin transport is ballistic. In the low frequency linear re- 
sponse limit, the time derivative of 3 S can be neglected. 
The spin current is thus constant in time, and the space 
dependence of the magnetization is a precession about 
the direction of the spin current [^|,^) . 

In the opposite limit of large temperature, the second 
term will dominate and the spin transport will be mainly 
diffusive. We retrieve the diffusion equation: 



3 



v 2 F d 2 M(X) + u 2 



k B T 
2i:h 3 'v 2 F 



d t M(X) = 0, 



(18) 



and we thus extract the spin diffusion coefficient: 
2Trh 3 vt 



u 2 k B T 



(19) 



This differs only by a factor of 8 from the value derived 
in§. 

The crossover temperature between diffusive and pre- 
cessional spin transport can be estimated as a function 
of the interaction strength and of magnetization 



knT, 



c roxsover 



M 



(20) 



IV. TIME EVOLUTION OF AN UNIFORM SPIN 
CURRENT 

Here we find the time evolution of a spatially uniform 
spin current with the only non-zero component on the z 
axis, injected into the system at t = 0. We do not use a 
perturbative approach as in the previous section, but we 
focus on the more general form in Eq.(^0|) and we try to 
study it both numerically and analytically. We assume I 



the initial magnetization is zero everywhere. This can be 
implemented by applying different chemical potentials on 
the right/left, spin up/down subspecies of particles. In 
particular we will choose chemical potentials /i for the 
spin up right movers and for the spin down left movers 
and respectively chemical potentials — /i for the spin down 
right movers and for the spin up left movers. The initial 
Fermi functions for the electrons are thus space indepen- 
dent and can be written as: 



= 0, a) 
fg(p,t = 0,x) 



flHp,t = o,x) = 

f l L l (j>,t = 0,x) 



1 





+ 1 


1 




e k B T 


+ 1 


1 






+ 1 


1 





e * + i 



f% ia2 (p, t = 0,x) = 0, for all a x ^ a 2 . 



(21) 

(22) 

(23) 

(24) 
(25) 



For this particular case, we can write simplified equa- 
tions involving only the diagonal components of the dis- 
tribution function matrix. From Eq.([Tc|) we derive 



± 



4:h VF L Jq Jq 

+ I fl 1 (Q)f L H L (p)-flHp) I ft l (Q)}, 



(26) 



The equations for the distributions functions of the 
left movers are obtained simply by interchanging L and 
R everywhere in the above formula. These equations can 
be solved numerically for finite temperatures when we 
impose the initial conditions from Eqs.d2l] - Eq). 



d t 3l = 



i k B T 
2irh 3 v 2 F ' 



(27) 



with the solution 



A. High temperature regime 

In the limit of temperatures that are much larger than 
the applied chemical potential, the numerical analysis in- 
dicates that the system evolves towards an equilibrium 
state described by the usual Fermi distributions. 

We can compare the numerical results with an analyt- 
ical estimate by noting that in this limit the quantum 
Boltzmann equations for the system are easily solvable. 
This is because at each momentum, the variations from 
the equilibrium distribution function are small, and we 
can describe the system by the equations derived in sec- 
tion III, specialized to the case of zero magnetization. 
From Eq.fJI) we derive: 



3l(t)=JU0)e- tV , (28) 



where V= 

This is confirmed by the numerical analysis. When we 
plot the logarithm of the current as a function of time 
we obtain straight lines, consistent with the assumption 
that the current is decaying exponentially. If we set all 
the constants to 1 the slope of the decay is just — u 2 T. 
We compare the numerical results with the theoretical 
estimates for an exponential decay (represented in the 
figure by the dashed lines) and the results are in very 
good agreement. 



4 




FIG. 1. The logarithm of the spin current plotted as a func- 
tion of time. The interaction strength has been taken equal 
to 0.3, the chemical potential [i = 1 and the curves denoted 
by 1, 2, 3 and 4 correspond to the temperatures of T = 4, 7, 
10 and 20 respectively. 



B. Zero temperature regime 



In the opposite limit of zero temperature (or fi ksT) 
the numerical analysis indicates strong deviations from 
the purely exponential decay. This can be understood 
as follows: at low temperatures the perturbations cannot 
be treated as small with respect to the equilibrium Fermi 
function. As an example, in FigJ| we plot the distribu- 
tion functions for the case we described in Eqs.(|2ll-p5|), 
at t = and our numerical predictions for the form of 
these distribution functions at later times. In the range 
\p\ < (x/hvF the difference from the equilibrium distribu- 
tions can be as high as 1, therefore making the pertur- 
bative approach of section III inconsistent. 

We note that instead of moving toward the unper- 
turbed Fermi distributions, in the region \p\ < fi/hvF 
the value of the distribution functions is a constant. At 
very large times this constant goes to 1/2. This type 
of behavior is entirely consistent with the physics of the 
system. In one dimension with linear dispersion, for the 
specific form of the interactions we consider, two interact- 
ing electrons can only transfer spin while the momenta of 
the electrons are preserved. Thus, the sum of the num- 
bers of spin up and spin down right moving electrons at 
each momentum has to remain constant in time. 
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FIG. 2. The distribution functions for a) spin up right 
movers, b) spin down right movers c) spin up left movers 
d) spin down left movers depicted at t — and at subsequent 
times. After sufficiently long times, the values for all the in- 
dividual distribution functions will be constant in the region 
—fi/hvF < p < /j,/hvF and equal to 1/2. Here we set all the 
constants to 1 and we chose /i = 2, T = 0.1 and u = 0.3 

Taking this observation into account, we can also solve 
analytically the equations of motion for the spin cur- 
rent in the limit of zero temperature. We define f R 

![/£_- /A 1 ] ^d/£ = /^ +/j 



Eq.(p6[) becomes 
dtf s R/L (p) = 

-fk/M 

+1[f S R/M 



Jj . In these new variables 



f & L/R(l)[fR/L(p)~llfR/L(p)f 



(29) 



/£(«)-„ l[ft/R(p)Y 

q 



/£/*(<?) -2 [[fl/R(q)} 2 fR /L (p)}, 



For the case defined by Eqs.(pT| -p5|), the form of the 
analytical solution can be inferred from the results of the 
numerical analysis. In Fig.|| we plot j°J { J L at t = and 
subsequent times as obtained from the numerical analy- 
sis. 

Using Eq. (p6|) , we note that dtj R j L = 0. Therefore 
these distribution functions are constant in time and at 
any time we can use for them the values we fixed at t = 0. 
For this particular case, 



/£(p) = 0(-p- 



Tivp ' 



G(-p- 



) 



and 



flip) 



9(p-rf-) 
Tlvf 



e( P + ^). 

nvp 



(30) 



(31) 
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Then Eqs.([29|) simplify, so that for \p\ < fi/hvpi 

u 2 r l 



u 2 (If 

dtfSR ' L ^ )= w^MlJ l ' R{q) 



+m /L (p)} 2 J fi/ R (i) - 2 j[f s L/R (q)?r R/L (p)}, 

while 

9tf R/L (p) = (33) 
for other values of p. 



J*(0) 



v / 2e*/ t ° - 1 ' 



(35) 



where J 2 (0) = 2^/2-kTivf and t = 2nh 3 Vp//j,u 2 . 

This result agrees very well with the numerical analy- 
sis for all values of the parameters in the system. Below 
we plot the numerical versus the analytical result, for a 
particular set of values for T, /i and u. 





FIG. 3. The distribution functions a)/£, b)/£, c)/£, d)/£ 
at t = and at subsequent times. The value for fj, has been 
set to 2 while the value for u is 0.3. We note that the "charge" 
distribution functions are time independent, while the "spin" 
distributions functions are constant on the interval (— /U, /*), 
and decrease to zero after long enough times. 

From Eqs.( |3^ , |33| ) as well as from the numerical anal- 
ysis plotted in Fig.|], we note that at t = 0, f R (p,t) + 
f[{p, t) = 0, which implies that f^{p, t) + f[(p, t) = at 
all times. Moreover we can pick a solution of the form 



Sr/lM 



±A(£),for \p\ < [i/hvF, 
0, otherwise. 



with A(0) = 1/2. Therefore the equation for A becomes: 



dX 

17 



ATrh 3 v 2 F 



(A + 4A 3 .) 



(34) 



Noting that 3 Z = J p [f s R (P, t) - f[{ P , t)] = 4\(t)ii/2irhv F , 
we obtain a time dependence for the spin current of the 
form: 




FIG. 4. The spin current is plotted as a function of time, 
for T = 0.1, fi = 2 and u = 0.3. Both the numerical and an- 
alytical result are given; the dashed lines corresponds to the 
estimate of Eq.(p5[), while the full line is the numerical result. 

We note that for short times it <C to) the form of the 



decay is a power law 3 z (t) 



_ J'(0) 



while for large 



enough times t 3> to we retrieve a purely exponential 
decay J 2 (t) = ^M e -^ 2t ° 



C. Intermediate temperature regime 



We can perform the numerical analysis for the inter- 
mediate regime of temperature when the temperature is 
comparable to the applied chemical potential. We study 
the crossover between the high temperature and low tem- 
perature regime by plotting the rate of change TZ of the 
log of the current at large enough time scales as a func- 
tion of ji or T. As seen in the previous sections, in the 
zero temperature limit we expect this to be a constant 
TZ = — l/2io ~ — fiu 2 /2, while in the high temperature 
limit 1Z = -V u 2 T 

In Fig.|5| a) we plot TZ as a function of T, while keeping 
the chemical potential fixed to fi — 1, and in b) we plot 
TZ as a function of of while keeping the temperature 
constant T = 1. 
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FIG. 5. The rate of change of the log of the spin current 
with respect to 

a) the applied chemical potential fj, when the temperature is 
held constant T = 1. 

b) the temperature T, when the chemical potential is held 
constant to fi = 1. 

The value for u is 0.3, 1/to = 0.09, and V = 0.09. 

As we showed before, when T is held fixed and [i is 
small R goes to a constant value V ~ —u 2 T = —0.09. 
Similarly, when /i is fixed and T is small R becomes - 



l/2t 



2 /i/2 = —0.045. The figures above also por- 



tray the crossover between the two regimes. 



V. CONCLUSION 

We analyzed the effects of backscattering interactions 
on spin transport in a one dimensional quantum wire. 
Using the quantum Boltzmann equation we wrote down 
the equations of motion for the spin current and mag- 
netization. We found that in the hydrodynamic limit of 
small perturbations, the backscattering processes gener- 
ate terms in the equation of motion responsible for bulk 
precession of magnetization as well as for spin diffusion. 
We computed the diffusion coefficient. 

We also analyzed both numerically and analytically the 
time evolution of a uniform spin current injected into the 
system. In the limit of high temperatures, this shows an 
exponential decay controlled by the values of the temper- 
ature and of the interaction strength. In the zero temper- 



ature limit, the analytic results indicate a more complex 
behavior characterized by a transition between a power 
law behavior at small times to an exponential decay be- 
havior controlled by the strength of the applied chemical 
potential and by the strength of the interactions, at large 
enough times. 
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APPENDIX A: 



1. Equations of motion for the imaginary time Green's function 

Given the non- interacting low energy effective Hamiltonian for a spin-i quantum wire described in Eq. ([!]): 

H = ^2 dx nv F [ipl a id x tp La - tp Ra id x ip Ra ], 

a J 

and the local interactions defined in Eqs.(||,||): 

H\ = -u f J dx.J R (x)J L (x) = -Uf^2 J dxip Ha (x)il) Ra (x)i()^(x)ip L p(x), 

a, P 

and 

H 2 = —u J dxJ R (x) ■ J L (x) = ~ ^ f dxif) Ra (x)a a pil) R p(x) ■ 2pl 1 (x)a 1 sipLs(x), 

we derive the transport equations for the spin. 
We start with the equations for tp ra : 

ihd t ^ ra {x,t) = [H(t),ipra(x,t)}, 

which implies 

ih[d t - Kv F (-l) er d x ]ip ra (x,t) = - 7;8 rR ^ ipj^jx, t)ip La (x, t)tp Rl (x, t) 

7 

+(j " u f )6 rR ^2ipl 7 (x,t)ip Ll (x,t)il; Ra (x,t) + L), 

7 

where r = R/L, e R ; = 1 and ex = 2. 

We define the imaginary time Green's functions to be: 



2{<*l ri :^2r 2 i*3r 3 ) ^r 4 . 



. <T[S^ ri M)^{X 2 )] > 
1 < T[S] > 

< T[S] > 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



where oi\,u 2 =T / -U r i — R/L, and Xi stands for Xi,ti. Also, T denotes imaginary time ordering, and S is given 



by 



S = exp 



-ip 



n EEr dtu%,{x)i,i r {x)i, ir ,(x) 



Using Eq. ( |A6| ) we can derive the equations of motions for the imaginary time Green's functions: 
%[d h - {-l)^v F d xl ]G^(X 1 ,X 2 ) + iJ2u^ a {X 1 )G^(X 1 ,X 2 ) = -iM{X x - X 2 )S aia2 S rir2 



(A7) 
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Because momentum conservation implies G^* 2 = 0, for simplicity we will denote G"^ a2 = G" lQ2 . 
We introduce the self energy for the right- and left- moving particles by the relation 

G-l ai a 2 ,y y \ ft — I&l0t2 ( V "V \ X^ a l a 2 { y V \ 
R/L \<*-li<*-2) — <^Q R / L \<*1,<*2) - ^ R / L l/VLj<"<"2j, 



(A8) 



where Go is the Green's function in the absence of the backscattering interactions between particles. 

The self energy can be computed perturbatively. We generalize the method of |l4| to include chirality and spin 
indices. In the collision approximation (second order perturbation theory in the interaction), the self energy is 



V£j2{X u X 2 ) = i^-G^(X 1 ,X 2 )S(X 1 - X 2 ) - ~(~ - u f )J2G a L J R (X ll X 2 )S(X 1 X 2 )5 a%c 



(A9) 



\ ' r«i n 2 I v y ^/-^«P ( y y \r<P<x (y y \ 
' ~^2 G L/R i^l > <*2 )G L/R \*2 , Xl)G R/L {Xl , X 2 ) 



a,f3 

u(u — 4u/) 
u{u — 4u/) 

16?i 2 



Er-mia / y y \/~iaf3 f y y \y^,f3a 2 j y y \ 
G L/i?^i' *2)G L/R {x 2 , x x )g r/l {a: 1 , x 2 ) 



a, (3 



E(~taia i y y \/~iaj3 iy y \W3a 2 / y y \ 
^RIL\^^2)^ l i r \X 2 ,Ax)(j l i r \Ax,X 2 ) 



E^a 1 a 2 /y y \/^a(3 iy y \(-if3a iy y \ 
'-'r/l \Xi,x 2 )Lj l/r [a 2 ,a:i)Lj l/r {a:i,a: 2 ). 



a, (3 



2. Equations of motion for real time Green's functions 

Following |l4| ] we perform the analytic continuation from the imaginary time to real time. We define real time 
Green's functions 



g? ia *(X u X 2 ) = -i< T[Vy ai (WL(* 2 )] >, 



(A10) 



a±ot2 > 



{X U X 2 ) =-i< ^1(^1)^(^2) >, 



We also define 



?? lQ2< (A' 1 , X 2 ) = l < $a a (Wr ai (*i) > 

gf(x,t)=gf(x,t;x,t), 



and we introduce the Fourier transforms: 



?^ </> (p, W ;x,t)= [ dx' f^'e^lTirtil^+^-l'-^) 



(All) 
(A12) 

(A13) 
(A14) 

We work in the collision approximation. We also assume that the space and time dependence of the Green functions 
of the system g(Xi,X 2 ) vary slowly with respect to the center of mass coordinates x = {x\ +x 2 )/2 and t = (ti+t 2 )/2 
|l4j . Thus the equations of motion for the real time Green's functions become: 

h[d t ± v F d x ]g%)f<(X) = ± £ || [g^ < (X)gr < (X) - g^ < {X)g R ^ < {X)] j 

±^ E / / l\g^ < ( Pl ,w 1 ;X)g R a2> (p 2 ,w 2 ;X)g^ < (p 3 ,w 3 ;X)g2 3ai> 
-(R <-> L)J - (<^>)^2tt5(pi +p 3 - p 2 - P4)2ttS(uji + uj 3 - uj 2 - 0J4) 
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,u(u — 4uf) \ 



16^ 



92 ia3< (Pi, w H X)gT ai> (P2, w 2 ; A , )^, 47< 3 , w 3 ; X)g] i a2> (p 4 , w 4 ; X) 
(R<~*L) - (<<->>) >27T(5(pi +P3 — P2 -P4)27r5(wi + cj 3 — w 2 - w 4 ). (A15) 



Here we denoted |f by J^, and similarly for ui. As before we chose X to denote (x, i). For simplicity we set 
the high energy cutoff in the momentum and frequency integrals at infinity, as the energy scale beyond which the 
Luttinger Liquid physics ceases to be valid is much larger than the other energy scales involved in the problem. 



3. The Quantum Boltzmann Equation 

We introduce the generalized Wigner distribution functions by: 



f^(p;x,t) = j dx , e- i P*'{-ig?^ < (x + j,t;x-^,t)}= j [g^<(p,u;;x,t)]. 

The equations for these functions can be derived from Eq.( A15| ) realizing that Jl4| |: 

g^fawiX) = f^(p;X)af(p,u;), 



g^ip^iX) = [5 a0 - f^(p;X)}af(p,u), 
a a L ^ /R (p, uj) = 2ttS(uj ± v F p). 



(A16) 

(A17) 
(A18) 
(A19) 



We thus obtain: 



(d t ± VF d x ) [ /s)2"(p)=±^iE / ifR 1 (pi)fr(P2)~fr i (pi)fR 2 (P2)} 

Jp Zn „ Jpi,p2 



± 



8h 2 VF 



E 



/f 7 (pi)[<W - f R a2 (P2)}fT a3 (P2)l^ - f2 3ai (pi 



7,q3,q 4 J Vl,Vl 

- /r i7 (pi)]/r 2 b2)[^ 4a3 - /« 4Q3 (P2)]/r Q4 (Pi) 



- (R <-> L) 



(A20) 



Using the same approximation of slowly varying distributions functions we can write these equations so they describe 
the change in the distribution functions at each momentum p: 



(d t +v F d x )f^(p) = -i[f^(p) J f2 a2 (q)~ J fZ 1J (q)fr(p)} 



(A21) 



8h z v F 



+ 1 fr(a)a a2 (p) 



fT 3 ip) / f a L 3ai {q)- / fV{q) 



2 / f^{q) 



r L ia3 {q)r R 3ai (p)-pR{p) 



ft^{q;X)-f^{p) 



- / f a R 1 \p)fr{q) 



f% 3a3 (p)~fZ 3a3 (q) 



{d t -v F d x )r L ia2 {p)^^\fr{p) 1 ir(q)- \ iT^r^)\ 



(A22) 



+ 



0*2 1 "JL 

an v F 



Vl ia2 {p) 



ir 3 (p) 1 fT a4 (<i)- 1 



2 / f£°°(q) 



fR ia3 {q)t L 3ai {p)-fl\p) 



+ 1 fsrwfT'tp) 



fR 3a3 (q)-fL 3a3 (p) 



fL 3a3 (p)-fR a3 (q) 



These are the generalized Boltzmann equations for a one dimensional spinfull electron gas. 
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